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Dyson hierarchical quantum ferromagnetic Ising chain 
with pure or random transverse fields 

Cecile Monthus 

Institut de Physique Theorique, Universite Paris Saclay, CNRS, CEA, 91191 Gif-sur-Yvette, France 

The Dyson hierarchical version of the quantum Ising chain with Long-Ranged power-law ferro¬ 
magnetic couplings J(r) oc and pure or random transverse fields is studied via real-space 

renormalization. For the pure case, the critical exponents are explicitly obtained as a function of 
the parameter a, and are compared with previous results of other approaches. For the random case, 
the RG rules are numerically applied and the critical behaviors are compared with previous Strong 
Disorder Renormalization results. 


I. INTRODUCTION 

The quantum Ising model is the basic model in the field of zero-temperature quantum phase transitions [I|. For 
the short-ranged case on hypercubic lattice in dimension d [I[, the dynamical exponent relating time and space via 
t = takes the simple value 


4T = 1 ( 1 ) 

This means that the time simply plays the role of an additional spatial dimension in the quantum-classical corre¬ 
spondence, so that the d-dimensional quantum Ising model is equivalent to the classical Ising model in D = d 1 
dimensions. It is also interesting to consider the long-ranged quantum Ising chain 

= (2) 

i i,j 


with uniform transverse field h and power-law ferromagnetic coupling 




d-J 


U-o- 


(3) 


in order to study how the critical properties depend upon the parameter cr > 0 0 . This quantum chain is equivalent to 
a 2D classical Ising model with nearest-neighbor coupling in the ’time’ direction and long-ranged power-law coupling 
in the ’spatial’ direction, so that the dynamical exponent z(a) is not given by the short-ranged isotropic value of Eq. 
[T] anymore, as a consequence of the strong anisotropy between time and space. In particular in the mean-field region 
0<(T<cr„ = 2/3 (see the reminder in Appendix]^, the dynamical exponent z, the correlation length exponent 
iz and the anomalous dimension 77 read Q 


< a < = 2/3) =1 

<cr <(Ju = 2/3) = i 

(0 < cr < CTu = 2/3) = 2 - cr 


(4) 


It turns out that the other problem of the dissipative short-ranged quantum spin chain is equivalent after integration 
over the bath degrees of freedom to a 2D classical model with nearest-neighbor coupling in the ’spatial’ direction and 
long-ranged power-law coupling in the ’time’ direction, where the parameter a now characterizes the bath spectral 
function J(a; —>■ 0) oc w®’ the value cr = 1 corresponds to the most studied Ohmic damping, whereas the region 

0 < cr < I corresponds to sub-Ohmic damping, and the region cr > 1 to super-Ohmic damping. The dissipative 
short-ranged chain is thus equivalent to the model of Eq. [5] after the interchange of space and time. So the dynamical 
exponent z', the correlation length exponent v', and the anomalous dimension 77' of the dissipative quantum chain 
are given hj z' = 1/z, v' = vz and z -|- 77 — 1 = ^ ~^ As a consequence, the values z' ~ 2 , tz' ~ 0.64 and 77' = 0 

measured via Monte-Carlo for the dissipative Ohmic chain [1, @ translate for the quantum chain of Eq. [5] into 


z£-^(a = l) =.0.5 

7zr-(a = l) =.1.28 

77£r(a = l) =.1. 


( 5 ) 
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Let us now consider the effects of randomness in the transverse-helds hi 


Tjrandom 

^LR 



(id) 


( 6 ) 


The relevance of a small disorder at pure quantum phase transitions needs to be discussed from two points of view 
0 : on one hand the Harris criterion Q or equivalently the Chayes et al inequality Q imply that the pure fixed 
point can be stable only if > 2/d = 2 here in the spatial dimension d = 1; on the other hand, the analysis of 

rare regions Q shows that it is important to compare the dimensionality djm of rare regions (here dun = 1 since 
the disorder is actually ’infinitely’ correlated along the time-direction) and the lower critical dimension = 1 

sufficient to obtain ferromagnetic ordering : the case dRR = corresponds to models where rare regions can play 

an essential role at criticality (see Q for more details). The random LR chain of Eq. [5] with a finite initial disorder 
has been studied recently via the Strong Disorder Renormalization (see for a review) : the main results 
(see also the related work |12j| in arbitrary dimension d) are the following critical dynamical exponent 

4DHr(^) = i+^ (7) 

and the essential singularity of the correlation length as the control parameter 6 approaches its critical value 9c 

( 8 ) 

corresponding formally to an infinite correlation length exponent 

(9) 

These properties should be contrasted with the short-ranged random Chain governed by an Infinite Disorder Fixed 
Point characterized by an infinite dynamical exponent 

4k = +00 (10) 

and two hnite correlation length exponents ^ and v'gRav™' — ^ [i3l' Note that the effects of disorder on 

the dissipative quantum chain mentioned above has been also much studied via Strong Disorder RG [3, [HI and via 
Monte-Carlo [31 j but here the problem is not equivalent to Eq. [SJ as a consequence of the columnar nature of the 
disorder along the time direction in each problem (the interchange of space and time discussed above for the pure 
case is not possible anymore). 

Since the obtained critical dynamical exponent of Eq. [3 is finite and not infinite as in the short-ranged case (Eq. 
[TOD, the Strong Disorder RG approach is not exact asymptotically, but only approximate. As a consequence, it seems 
useful to analyze the critical properties with another approach in order to compare the results. For the random 
Short-Ranged chain, the seld-dual block RG procedure first introduced for the pure chain [3 has been found recently 
to be able to reproduce the Fisher Infinite-Disorder fixed point (l8l - [20| : this shows that the same block RG rules 
can lead to conventional critical behavior for the pure chain and to Infinite disorder critical behavior for the random 
chain depending on the initial condition of the RG flow. The aim of the present paper is to study via some block 
renormalization the Dyson hierarchical analog of the pure and the random long-ranged chain of Eq. [2] and Eq. [6l 
The paper is organized as follows. In section im we introduce the Dyson hierarchical quantum Ising model and 
derive the renormalization rules. In section uni the RG equations for the Dyson analog of the pure long-ranged chain 
of Eq. [3 are solved analytically. In section IIVI the RG rules for the Dyson analog of the random long-ranged chain 
of Eq. [6] are studied numerically to obtain the critical properties. Our conclusions are summarized in section El 
Appendix contains a reminder on the mean-held theory for the pure long-ranged chain. 


II. RENORMALIZATION RULES FOR THE DYSON HIERARCHICAL MODEL 

A. Dyson hierarchical version of the Long-Ranged Quantum Ising Chain 

In the held of long ranged models, it is very useful to consider their Dyson hierarchical analogs, where real space 
renormalization procedures are usually easier to dehne and to solve as a consequence of the hierarchical structure. 
The Dyson hierarchical classical ferromagnetic Ising model [2l| has been much studied by both mathematicians (3 - 
[H| and physicists [Hl - [30| . More recently, Dyson hierarchical versions have been considered for various disordered 
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and spin-glasses 


or quantum like Anderson 


systems, either classical like random fields Ising models 
localization models 

Here we introduce the Dyson hierarchical analog of the Long-Ranged Quantum Ising Chain of Eq. [6] as follows. 
The Hamiltonian for 2" quantum spins can be decomposed as a sum over the generations fc = 0,l,..,n—1 


n — 1 




(k) 

( 1 . 2 ” 


k=0 


The Hamiltonian of generation k = 0 contains the transverse fields hi and the lowest order couplings 

2” 2"“1 

<2°) - E 


2 = 1 


2=1 


The Hamiltonian of generation k = \ reads 


'^(1,2^) — ^ ^ ^ [a^ 42—3^42—3 M42—2'^4i_2] [/^42—l<^4i_i + f^4.iO'4i\ 


2 = 1 

the Hamiltonian of generation k = 2 reads 

271-3 


Hi 


(E") “ “ E [^k8^-70^si-7 + MSi-eO'si-e + hSi-50'si-5 + M8i-40’8i-4] 

i=l 


( 11 ) 


( 12 ) 


(13) 


(14) 


X [^8i-30'gi_3 -I- fi8i-2<^8i-2 + MSi-lCgi-i + hSiC^Si] 
and so on up to the last generation k = n — 1 that couples the two halves of the system 


(„-i) _ _ 
"(1,2") “ 


277-1 

277 

E 

E 

2=1 

_i=2"-i-|-l 


(15) 


The transverse helds > 0 in Eq. [T^can be either uniform or random. The magnetic moments of the spins are 
set initially to unity 


At* = 1 


(16) 


but we have introduced them because they will be generated by the renormalization procedure described below. 

The > 0 are given ferromagnetic couplings as a function of the generation k. To mimic the power-law behavior 
with respect to the distance r of Eq. [3] 


J(r) = 


-.1 + (T 


we consider the following exponential behavior with respect to the generation k 


= 


1 


(2fc) 


l+rr 


_ 2“(n-<’’)fe 


(17) 


(18) 


At the classical level, the energy cost with respect to the ground state of a Domain-Wall between the first half-system 
having 5'“ = 1 for 1 < A < and the second half-system having S'® = — 1 for ^ -|- 1 < j < L scales as 


E^^{L) oc 


(19) 


In the region ct > 1 where it decays with the system-size L, it is clearly different from the Long-Ranged model of Eq. 
[2] that at least contains a constant term coming from the nearest-neighbor coupling between the two halves. On the 
contrary in the region cr < 1 where Eq. [19] grows with the system-size L, one may expect that the Dyson hierarchical 
version is an appropriate approximation of the Long-Ranged model. The case cr = 1 is at the border line, since the 
Domain-Wall cost remains constant for the Dyson model (Eq. [T31) . whereas it grows logarithmically in L for the 
Long-Ranged model. In the following, we will thus focus on the interval 

0 < cr < 1 (20) 

In particular we should stress that whereas the limit a —>■ -l-oo of the Long-Ranged model of Eqs |2| and |6| yields the 
corresponding Short-Ranged models, the limitation cr < 1 for the Dyson hierarchical versions does not allow to recover 
the Short-Ranged models in this limit. 
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B. Diagonalization of the lowest generation Hamiltonian 


The Hamiltonian of generation k 

H(2i-l,2i) = 


= 0 of Eq. [12] is the sum of the independent two-spin Hamiltonians 
~^2i—lO'2i—1 ~ h2i<J2i ~ k‘2i—lt^2i(^2i—1^2i 


( 21 ) 


1. Diagonalization in the symmetric sector 


Within the symmetric sector, the diagonalization of the Hamiltonian of Eq. j^in the basis 
^(2i-l,2z)l + + > = ~(^2z-l + ^2i)| "b "b > — lJ-2i-lk‘2i\ -> 

H(2i-l,2i) I-> = lJ‘2i-llJ‘2i\ "b “b > —(/l2z-l “b /l2i)|-> 

leads to the two eigenvalues 


Afj = (J(°V2i-lM2i)^ + {h2i-l + h2iY 

= +\J + {h2i-l -b h2iY 

with the corresponding eigenvectors 

|Af^" > = cos6»|j -b -b > +sin6»fj - - > 

|Af+> =-sin6»|j-b-b >-bcosdfil - - > 


in terms of the angle 9$ satisfying 


cos(ei^\ 


\ 


sin(dD) = 


\ 


1-b 

h2i-i+h2i 

yj fl2i-lfJ'2i)‘^ + {h2i-l-\-h2i)^ 

2 

X 

h2i-l-\-h2i 




( 22 ) 


(23) 


(24) 


(25) 


2. Diagonalization in the antisymmetric sector 


Within the antisymmetric sector 

^{2i-l,2i) \ -^ “ ~(^2z-l ~ h2i)\ H-> — l^2i-lt<'2i\ -b > 

^{2i-l,2i)\ -^ ^ “ —J^^^tJ-2i-ltJ-2i\ H-> +(^2z-l ~ ^2i)|-b > 

the two eigenvalues read 

^2i = ~\/92i-lk-2iY + (^2i-l — ^2i)^ 

^2^ — ~^\J 92i-lk‘2iY + (^2i-l — ^2i)^ 

with the corresponding eigenvectors 

|A^“ > = cos 0^,1-b->-b sin --b > 

|A^j+ > = -sind^jl -b - > -bcos0^j| - -b > 


in terms of the angle 9a satisfying 


cos(0^i) 

sin(0^i) 


1-b 

—/J'22 


2 

1 

h2i-l—h2i 



2 


\ 




(26) 


(27) 


(28) 


(29) 
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C. Introduction of the renormalized spins (jR(2i) 

For each two-spin Hamiltonian of Eq. [211 wish to keep the two lowest states among the four eigenstates 

discussed above, and to label them as the two states of some renormalized spin <jR(^2i) 

\'^R(2i) = + > = > 

\'^R(2i) — ~ > = \^2i > (30) 

It is convenient to introduce the corresponding projector 

^2i = kl;(2i) = + X WR(2i) = +1 + Wr{2i) = ~ X ^R(2i) = “I (31) 

as well as the spin operators 

^Ri2i) = Wr{2i) = + X Wr{2i) = +1 “ klj(2i) = “ >< ®'fl(2j) = “I 

^fi:(2i) = Wr(2i) = + X Wr(2i) = “I + klj(2i) = “ X 0-^(2*) = +1 (32) 


D. Renormalization rule for the transverse fields h 


R(2i) 


The projection of the Hamiltonian of Eq. |2T]is given by 

^2i^{2i-l,2i)P2i ~ ^2i \^2i ^^^2i I + -^21 l'^2i ^^^2i 


1^27 


= e.R{2i)P2i - hR(2i)(^R(2i 


\S- _ 

^2i -^21 


^R(2i) 


where the renormalized transverse fields read 
hli{2i) = 


— ^2i ^2i 


_ 'sj (•/^°V2i-lAi2i)^ + (^2i-l + ^2i)^ — \/(J(°V2i-lM2i)^ + (^2i-l — /i2i)^ 

“ 2 
_ 2/l2i_l/l2i 

\l (j('’V2i-l/^2i)^ + (^2i-l + ^2i)^ + \/ (J^°V2i-lM2i)^ + (^2i-l — /i2i)^ 

and where the contribution to the ground-state energy of this projection reads 

— ^2i + ^2i 

^R(2i) = - 7^ - 


yj ^^.2i-l^J■2iY + (/i2i-l + ^2i)^ + \/( J^°V2i-lM2i)^ + (^2i-l — /i2z)^ 


(33) 


(34) 


(35) 


E. Renormalization rule for the magnetic moments Rii(2i) 


The projections of the operators 

P 2 i^ 2 ^-lP 2 ^ = [sin( 6 'fj cos( 6 »^i) -f cos( 6 »fj sin( 6 »^J] cr^( 2 i) 

P 2 r^ 2 iP 2 i = [cos(0|i) cos( 6 »^J -f sin( 6 »|j sin( 6 »^J] ( 7 ^( 3 ^) 

yield the following renormalization rule for the magnetic moment 

^^R( 2 i) = [sin( 6 »fi) cos( 6 »^,) -f cos( 6 »|j sin( 6 »^J] ^ 2 z-i + [cos( 6 »|,) cos( 6 »^,) -f sin( 6 »D) sin( 6 »^J] ^MR( 2 ^) 


(36) 


— ^J‘2i-l 


\ 


1-f 


y/ fl2i-lfi.2iY + {h2i-l+h2i)'^\/ ( ^21-1 M2i - 1 —/l2i ) ^ 


+ M2i^ 


1-f 




2 


(37) 
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F. Renormalization of the Hamiltonians of generation fc > 1 

The Hamiltonian of generation k = 1 oi Eq. [13] is projected onto 
-1 \ /2"-i \ 2"-^ 

n ) ^(1 2'*) ( IT 1 = “ ^ Mfi(4i-2)0’^(4i-2) Cr^(4i) 

i=i / V*=i / 

The Hamiltonian of generation fc = 2 of Eq. [151 is projected onto 


(38) 


IT ^2i ^(1,2") IT ^^21 - “ X! Mfl:(8i-6)0'fl(8i-6) + Mfl(8j-4)0'ii(8i-4) Mfl(8j-2)0'fl(8i-2) ++Mfl:(8i)0'^(8i)(39) 


and so on up to the Hamiltonian of last generation fc = n — 1 of Eq. [T5] that is projected onto 


n ^27 n ^27 = 




)^R{2i) 


X! ^«(2i)': 


fl(21) 


j=2"-2 + l 


(40) 


In conclusion, once the renormalization of the magnetic moments has been taken into account, the couplings 
are just translated by one generation for fc > 0 


jR(k) ^ j(fe+l) 


(41) 


III. DYSON HIERARCHICAL CHAIN WITH UNIFORM TRANSVERSE FIELD 

The pure Dyson Quantum ferromagnetic Ising model corresponds to the case where all transverse fields hi coincide 

h^ = h (42) 

The mean-field theory for the Long-Ranged model of Eq. [Hi (see the reminder in Appendix [^ is expected to be 
also valid for the Dyson hierarchical version in the same region 0 < cr < CTu = 2/3. The main goal of the real-space 
procedure described below is to study the critical properties in the non-mean-field region cju > 2/3, but since a is just 
a continuous parameter in the RG rules, we will also mention the results for 0 < cr < CTu = 2/3. 


A. Pure renormalization rules 


The renormalization rules derived in the previous section simplify as follows 
(i) the renormalized transverse field of Eq. 1341 reads 


h^ = 


2h^ 


J(o)/ 42 -t v'(J(°V^)^ +4fc2 

(ii) the renormalized magnetic moment of Eq. [37] reads 


(43) 


hR 


— 1 -|- 


J(0)^2 


v7J(oVF+^ 


(44) 


B. RG rules deep in the paramagnetic phase h —>■ -too 

Deep in the paramagnetic phase h —>■ -too, the transverse field remains unchanged upon RG (Eq. 

h 

h—¥-\-oo 


(45) 
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whereas the magnetic moment evolves according to (Eq. 1441) 

~ (46) 

h—¥-\-oo 

For a length L = 2" obtained after n RG steps, the magnetic moment reads 

/r(L = 2") ~ (\/2)" = yi (47) 

h—¥-\-(X) 

As a consequence, the effective ferromagnetic coupling between two such magnetic moments scaling as 

= 2")m(A = 2”) ~ = (48) 

h—¥-\-oo 

becomes smaller and smaller with respect to the transverse field of Eq. 1451 so that the paramagnetic fixed point is 
attractive for all cr > 0. 


C. RG rules deep in the ferromagnetic phase h ^ 0 


Deep in the ferromagnetic phase ft. —>■ 0, the magnetic moment evolves according to (Eq. 

~ 2^ 
h—¥0 

i.e. the magnetic moment grows linearly with the length L = 2" obtained after n RG steps 

ii(L = 2’") ~ 2" = L 

h^O 


(49) 


(50) 


This corresponds as it should to the maximal magnetization per spin m = = 1. As a consequence, the effective 

ferromagnetic coupling between two such magnetic moments scales as 


= 2”)Ai(T = 2”) 




in agreement with the classical Domain-Wall energy of Eq. [TO] 

Deep in the ferromagnetic phase ft —>■ 0, the RG rule for the transverse field becomes (Eq. 


ft« 


ft2 


ft2 


h -,0 J ( 0)„2 j { 0 ) 

After iteration over n RG steps corresponding to the length L = 2" 

1 , 2 " 

TIN _ . fi" 


h{L = 2”) = ft" ~ 


ft^” 


/i->0 7(’^-3 )'i4 7 t(0) \2"-i TT"-! / xlfe) \2 "-i- 

We// ) We// ) "•We//! llfc=oWe//l 

or equivalently in log-variables using = 2fc(i-<^) (Eq. [51]) 

n—1 

In h{L = 2") 2" In ft - ^ 2"-^-'= ln( J^j}) 


= 2 " 


fc=0 
n —1 


lnft-^2-l-'=fc(l-cr) In2 


fc =0 


= 2" [Inft- (1 - cr) ln2 (l - (n-h 1)2“")] 

= L [In ft — (1 — cr) In 2] -|- (1 — cr) (In L -|- In 2) 


At leading order, the renormalized transverse field is thus attracted exponentially in L towards zero 

L 


h{L) 


ft-rO V 2^ 


(2L) 


1 —(T 


(51) 


(52) 


(53) 


(54) 


(55) 


so that it becomes smaller and smaller with respect to the effective ferromagnetic coupling of Eq. |5TJ i.e. the 
ferromagnetic fixed point ft = 0 is attractive. 









D. RG rule for the control parameter K 


It is clear that the important parameter of the model is the ratio between the effective ferromagnetic coupling 
and the transverse field h 


K = 


h 


After one RG step, this control parameter is renormalized into (Eqs [43] and [44] ) 


(56) 




hR ./IOT4 


= MK) 


(57) 


The attractive paramagnetic fixed point corresponds to AT = 0, whereas the attractive ferromagnetic fixed point 
corresponds to K = +oo. From now on, we focus on the unstable fixed point between them and on its critical 
properties. 


E. Location of the critical point A'c((j) 

The critical point Kc corresponds to the non-trivial fixed point Kc = iftaiKc) of the RG rule of Eq. [57] leading to 

{2"^ - K,)./Kfn = Kl + 2 (58) 

So keeping the condition Kc < 2'^, we may take the square to obtain the cubic equation 

Kl - 2^-^Kl + AKc - 2(2*^ - 2"'") = 0 (59) 

In terms of the standard Cardano notations for cubic equations 

qcr-i 


and 


p = 4 — 
q 


q2 p3 


27 


= 2^"'^ - 
= (2-'^ + 


2 ' 7‘+2 2 ^' 


<7-2 


27 




u+ = 


q2 p3 

T ^ 


= 1^-2-'^- 
u_ = -sgn(p) 
= -sgn(2 - 


2<7+i 


-)3(T—3 


27 


+ (2- + 2'^-i)^/l + - 


q 

/g2 



V 4 

^ 27 


In 3 
2 In 2 


-a) 


-2- 


2cr+l 


-3 


27 


-(2-‘^ + 2'^-^)^/l + - 


one obtains that the only real solution of Eq. [55] reads 


(60) 


(61) 


Kc{^) = —5-h u+ + u_ 


(62) 
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It is a growing function of cr with the values 


K,{a ^ 0) 

= CT In 2 + 0(1 

Kc{<y = i) 

= 0.177439 

KcA = i) 

= 0.364946 

KcA = ^) 

= 0.497043 

K,{a = I) 

= 0.783243 


(63) 


F. Correlation length exponent v{o-) 


The correlation length exponent v measures the instability of the linearized RG flow around the critical point. It 
is determined by the derivative of the RG flow of Eq. [571 


2 


^'AKc) 


2{2 + K,s/W+A 

(if2+4) 


(64) 


This yields that ^{cr) is a decaying function of tr with the values 


v{a —^ 0) 


:^- + 0(l) 
a 



= 4.4406 


-A = ^) 

= 2.45I3I 



= 1.9602 


v{a = I) 

= 1.482 

(65) 


G. Magnetic exponent x{a) 

At criticality K(.^ the RG rule for the renormalized magnetic moment of Eq. | 




— /iV^i /1 “f 




/f^F+4 


yields after n RG steps corresponding to the length L = 2" 


;,(L = 2”) V2J1 + 


Kr 




= L 


1 — X 


with the magnetic exponent 


x{a) = 1 — 




In 2 

1 — cr ln(Ar^ + 4) 


( 66 ) 


(67) 


4 


81n2 


( 68 ) 
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It is a decaying function of tr with the values 



1 a 

x{a —>■ 0) 

" 2 " 4 ^ 

= 1 ) 

= 0.438914 

= i) 

= 0.380907 

2, 


^('"= 3 ) 

= 0.344141 

x{a = 1) 

= 0.275732 


O(a^) 


(69) 


H. Dynamical exponent z(a) 


At criticality Kc, the RG rule for the transverse field h of Eq. 23] 

h^ = h - ^- 

K. + ^/Kf+l 


(70) 


yields after n RG steps corresponding to the length L = 2" 


h{L = 2”) = h 


with the dynamical exponent 


In 


z{a) = 




= hL 


(71) 


) (7 — 1 + 4) 


In 2 


41n2 


(72) 


By construction at criticality, the dynamical exponent z also describes the scaling of the effective ferromagnetic 
coupling 


L-^ = JeffiL) = L- 


(73) 

i.e. it is directly related to the magnetic exponent x of Eq. |33| 


z[a) = a - 

-1-1- 2a;(cr) 

(74) 

with the values 



z{a —?> 0) 

= f + o(.^) 



= 0.127828 


z{a=\) 

= 0.261814 



= 0.354949 


z(a = 1) 

= 0.551463 

(75) 


I. Correlation exponent ? 7 (cr) 


Within the real-space RG procedure, the two-point spatial correlation function scales as the square of the magne¬ 
tization at the scale where the two spins are merged into a single renormalized spin 


C{L = 2’^) =< alal >oc 



oc L 


-2x 


(76) 
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As a consequence, the exponent ry of the correlation function defined by 


Civ) 


^-(d-2+z+r) 


after taking into account Eq. O takes the simple explicit value 

rj{a) = 2 — a 


(77) 

(78) 


that coincides with the mean-field value of Eq. IA211 The fact that rj keeps its mean-field value (2 — cr) even in the 
non-mean-field region until it reaches its Short-Ranged value rjsR is expected for the Long-Ranged quantum Ising 
model 0 for the same reasons as in the classical case 


J. Ground state energy 


In the pure case, the renormalization rule of Eq. |35]for the contribution of a single block to the ground state energy 
simplifies into 


gr = 


-p y'(j(0)^2)2 +4/^2 


(79) 


Eor a chain of L = 2" spins, the ground state energy can be obtained by summing all the contributions of the 
successive projections for the blocks up to the last step where the single remaining spin has for contribution {—hRn) 


E^^{L = 2”, K) = --^gr - —eR2 - —e/ja - ... - e^n-i - hRr, 


n—1 


= -E 


L 


fc=0 


21+fe 


e_Rfe - hRn 


(80) 


In the paramagnetic limit J = 0 = K, where the transverse fields are not renormalized (Eq. |45|) . one recovers as it 
should the contribution {—h) for each independent spin 

n—1 


E^,Ul = 2\J = 0) =-Y^^h-h = -hL 


k^O 


(81) 


In the ferromagnetic limit h = 0,K = -|-oo, where the effective ferromagnetic coupling between two such magnetic 
moments scales as EqlHIl one obtains 

n — 1 


E 


,GS 

ferro v 


(L = 2"> = 0) = - ^ = - ‘ (L - £>-') 


1 


k=0 


2(1-2-'") 


(82) 


as it should : the coefficient of the extensive term corresponds to the half of the sum of the couplings starting from 
one site, for instance the first one, i.e. (X^iS whereas the correction to extensivity in represents the 

’missing’ couplings of the region j > L. 

At criticality K = Kc-, where the renormalized transverse-fields involve the dynamical exponent z (Eq. [7T]and Eq. 
EH), the ground state energy reads 

E2f,.(L = 2". A-J = - ”f Ail±LpiE±3 _ k,. 




n—1 


= -^E 

= —h 


21+fe 
L 


- h2- 


2l + k 

k=0 
22z 


2-zk2z _ ^2- 


(2^ - 1)^ , 
L + Xr-- -L-^ 


■'21+^ - ' 21+"’ - 1 
i.e. the finite-size correction to the ground state energy per spin 

EgUL = 2^,K,) 


= —h 




21 +"= - 1 21 +" - 1 


(83) 

(84) 


is of order L ^ ^ as expected. 
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K. Discussion 

In summary, the real-space RG procedure for the pure Dyson hierarchical quantum Ising model has the advantage 
of being explicitly solvable for the various observables and to yield reasonable approximated values for the critical 
exponents. Besides the simple anomalous dimension rj = 2 — a (Eq. 1781) reproduced both in the mean-field region 
and in the non-mean-field region as it should, it is interesting to consider some specific values of a to compare with 
previous approaches : 

(i) for the non-mean-field value cr = 1, the dynamical exponent z{a = 1) ~ 0.55, and the anomalous dimension 
r]{a) — 1 are close to the Monte-Carlo measure quoted in Eq. [5] of the Introduction, even if the correlation length 
i/(cr = 1) = 1.48 is less good with respect to Eq. [5j 

(ii) at the upper critical value au = 2/3, the dynamical exponent z{a = |) ~ 0.35 and the magnetic exponent 
x{a = |) ~ 0.34 are close to the expected values of Eos IA25] and IA261 even if the value for the correlation length 
exponent = |) ~ 1.96 is again not so good with respect to Eq. IA25I 

(hi) at first order in the expansion near cr —>■ 0, the critical exponents ir(cr), x{a), z{(t) and the critical point location 
Kc{(j) found above by the real-space RG procedure actually coincide with the mean-field values recalled in Appendix 
El Eurther work is needed to understand if the renormalization procedure can be modihed to describe correctly the 
whole mean-field region 0 < cr < CTu = 2/3 and in particular its anomalous finite-size scaling properties. Indeed, 
the anomalous finite-size scaling properties of the mean-field region of classical statistical physics models has a long 
history that has been re-interpreted recently (see (dsl - l^ and references therein), and it would be very interesting to 
see how it can emerge explicitly within some real-space RG. Eor the present case, the mean-field critical properties 
can actually be reproduced by the real-space RG flow for the control parameter K = (instead of Eq. [ST]) 


K 


R 

MF 


= 2 " 


K 


1-K 


(85) 


with the critical point location Kc{<j) = 1 — 2 ^ and the correlation length exponent v{a) = More precisely, the 
corresponding rules for the transverse field and the magnetic moment read (instead of Eqs |43] and 1441) 


h^p = hVT^ 


^^MF — 


y/T^ 


( 86 ) 


corresponding to the dynamical exponent z{a) = § and to the magnetic exponent x(a) = 5 — f • Within the real-space 
RG perspective, it is not clear to us what arguments should be used to justify that this mean-field RG flow has to be 
preferred to describe the critical point in the whole region 0 < cr < cr^ = 2/3. 

IV. DYSON HIERARCHICAL CHAIN WITH RANDOM TRANSVERSE FIELDS 


In this section, we consider the Dyson hierarchical chain of Eq. [TT] where the transverse helds hi are independent 
random variables drawn with the flat distribution over the interval [0, W] 


P{h) = 


e{o<h< w) 
w 


(87) 


which is the distribution previously used in the Strong Disorder RG study 0 . The parameter W which represents 
the typical scale of the initial transverse fields is thus the control parameter of the transition. It will be actually more 
convenient to use 


0 = In W 

as in the previous Strong Disorder RG study 0- 


( 88 ) 


A. RG rules deep in the ferromagnetic phase W ^ 0 

When the transverse fields {h 2 i-i,h 2 i) are small with respect to the ferromagnetic term the RG rule 

for the magnetic moment (Eq. I37p is simply additive 

hR(2i) = M2i-1 + 1^21 (89) 
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As a consequence, after n RG steps corresponding to the length L = 2”, the magnetic moment is given by the number 
of spins 


KL) = L 


(90) 


and the effective ferromagnetic coupling between two such magnetic moments scales as 

jW = = 2")/r(L = 2’") ~ 

in agreement with the classical Domain-Wall energy of Eq. [121 

The renormalization rule for the transverse field (Eq. 13411 becomes 


hR(2i) 


h2i-lh2i 


h2i-ih 

rCO) 

^eff 


2i 


(91) 


(92) 


Upon iteration in log-variables, one can follow the same steps as in the pure case fEos 1551 and [Mil and one obtains 


L 

Inh(L) ~ In hi — L(1 — cr) In 2-|-(1 — tr) (In L-|-In 2) (93) 

i=l 

For large L, the Central Limit theorem yields that the typical value 

\nhtyp{L) =lnh(L) 

~ L (In hi — (1 — a) In 2) -|- (1 — cr) (In L -|- In 2) (94) 

becomes smaller and smaller with respect to the effective ferromagnetic coupling of Eq. 1911 i.e. the ferromagnetic 
fixed point is attractive. In this ferromagnetic phase, the variance grows linearly (Eq. I93|) 

Uar[lnh(L)] = LUar[lnhi] (95) 


B. RG rules deep in the paramagnetic phase W -too 


When the ferromagnetic term J^°V 2 i-ih 2 i is small with respect to the transverse fields (h 2 i_i, h 2 i) the RG rule for 
the transverse field (Eq. [34]) becomes 


2h2i ih2i 


i .'• /7 7 \ 

R(2^) ^ T- —r -TTT-TT = mm(h 2 i_i,h 2 U 

h2i-l + h2i + |h2i-l — h2i| 


(96) 


and the RG rule for the magnetic moment (Eq. |37|) reads 


hfi(2i) 


1 - sgn{h 2 i-i - h 2 i) 


h2i-l + 


1 -t sgn{h2i-i - h2i) 


-l^2i 


— S{h2i-1 < h2i)g2i-l + S{h2i-1 > h2i)g.2i 


(97) 


In this paramagnetic limit, the renormalized spin is thus simply the spin with the smallest transverse field that survives 
with its properties, whereas the spin with the highest transverse field is decimated and disappears. 

So after n RG steps corresponding to the length L — 2^, the magnetic moment has simply kept its initial value 

= 2") =1 (98) 


so that the effective magnetic coupling also keeps its bare value 

Jeff{L) = J^P'> = 


(99) 


The renormalized transverse field is given by the minimum value among the L = 2" initial random variables 

/i(L = 2”) = min(/ii, ^ 2 ,.., hi) (100) 
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Its probability distribution PL{h) reads thus in terms of the initial distribution Pi{h) of Eq. [53 


PL{h) =LPiih) 


n-\-oo 


1 L-1 


dxPi {x) 


= L 


9{0<h< W) 
W 


W -h 
W 


L-l 


( 101 ) 


For large L, it converges towards the exponential distribution (as a particular case of the Frechet distribution) 


J^L(h) 




L—>-+oo IF 


The characteristic scale 


W 

hyp(P) = -J- 


( 102 ) 


(103) 


decays, but remains bigger that the effective ferromagnetic coupling of Eq. 11001 so that the paramagnetic fixed point 
is attractive. 

Th e ap proximate renormalization rules derived here in the limit W —> +oo are similar to the Strong Disorder RG 
rules |10l| . However the numerical study described below yields that the paramagnetic phase near the critical point is 
described by other scaling behavior. 


C. Numerical study of the Critical Point 

To study the critical properties, we have used two methods : 

(i) on one hand, we have generated Ug = 4.10^ disordered samples of the Dyson chain containing N = 2^"^ ~ 1.68x 10^ 
spins, corresponding to n = 24 generations, and we have applied numerically the RG rules of Eos [Ml and 1571 

(ii) on the other hand, we have used the so-called standard ’pool method’ : the idea is to represent the joint 

probability distribution P„(/i, /i) of the renormalized transverse field h and of the magnetic moment /r of a renormalized 
quantum spin at generation n by a pool of M realizations where i = 1,2,...,M. The pool at generation 

(n -I- 1) is then constructed as follows : each new realization is obtained by choosing 2 ancestors spins 

at random from the pool of generation n and by applying the renormalization rules of Eqs [Mj and 1371 The numerical 
results presented below have been obtained with a pool of size M = 3.10^ iterated up to n = 100 generations. The 
pool method allows to study much bigger sizes and statistics, but introduces some truncation related to the size of the 
pool in the probability distributions. We have checked that in the critical region, the results of (ii) are in agreement 
with (i) for the common sizes. 

In both methods, for each renormalization step corresponding to the lengths L = 2"^, we have analyzed the statistical 
properties of the renormalized transverse fields and of the renormalized magnetic moments fiL. More precisely, we 
have measured the RG flows of the corresponding typical values 

In = In hi 

In^^^ =ln/ri (104) 

as a function of the length L for various values of the control parameter 0 of the initial distribution of Eq. [53 We 
have studied the critical properties for the two values cr = 2/3 and a = 1 with the following critical point 

0^ (^<7 = - 2.80725 

0c (cr = 1) ~ 2.14875 


(105) 
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Figure 1: Case o' = | 

(a) RG flow of the magnetization rn^ = —^ in log-log plot for various control parameter 6 

(i) in the ferromagnetic phase 6 < 6c, the flow converges towards a finite magnetization rriacid) ; 

(ii) in the paramagnetic phase 6 > 6c, the slope is (—1/2) ; 

(iii) at criticality, the slope = 2/3)) ~ —0.5 is not distinguishable from the paramagnetic slope. 

(b) ln(mcx3(0)) as a function of ln(6c — 6) \ the slope corresponds to the exponent 13{a = |) ~ 3 


D. RG flows of the magnetic moment jj,L 

. tup 

On Fig. [11(a), we show the RG flows of the magnetization (Eq. 110411 as a function of the length L for 

various control parameter 9. They can be summarized as : 


mT\o<ec 

— '^cso(^) 



L—¥-\-CiO 



oc 



L — ^~|“00 


mT\B>6c 

oc A{9)L-i 

L—>+oo 

(106) 


On the ferromagnetic side, we obtain that the asymptotic finite magnetization vanishes as a power-law (see Fig. [T] 

(b)) 

moo(0) oc (0c(107) 
and we measure the same value for a = 2/3 and cr = 1 

= ^) - 3 ^(cr = 1) (108) 

On the paramagnetic side, it is not clear that the amplitude A{6) diverges, and correspondingly at criticality, the 
magnetic exponent x(jj) cannot be distinguished from the paramagnetic phase value 1/2 

a;(cr = -) ~ 0.5 ~ x{a = 1) (109) 

As a consequence, the finite-size scaling matching between Eq. 11081 and 11091 corresponds to the finite-size correlation 
length exponent vps = § 

2 

’^Fs{o- = -) ~ 6 ~ PFs{o- = 1 ) 


( 110 ) 
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E. RG flow of the effective coupling Jl 

In the present real-space RG framework, the effective coupling is just the ’slave’ of the magnetic moments, and 
the typical value is given by 


jtyp ^ ^ 2 , 1 —( 111 ) 

The translation of Eq. 11061 corresponds to the behaviors 

JL^\e<e, cx 

L—^ + CXD 

T^yPl rl — (T— 2 x r—z(( 7 ) 

Jl \e=e, « = L ^ > 

L —^~|“Oo 

Jm9>9. , « L-'^A^iO) ( 112 ) 

L^+oo 

Eq [109] yields that the critical dynamical exponent is given by 

z{a) = tT—l-|-2x~cr (113) 

F. RG flow of the renormalized transverse fields 




Figure 2: RG flow of the transverse fields for cr = | 

(a) ln(-ln/i))'*’) as a function of In L : 

in the ferromagnetic phase 9 < 6c, the slope unity allows to measure the correlation length ^h{9) of Eq. 11151 

(b) ln^h(0) and In^^idthiS) as a function of \n{6c — 9) : the slope corresponds to the value i^((j = |) ~ 6. 

At criticality, the typical renormalized transverse field involves the same dynamical exponent as in Eq. 11131 

oc (114) 

L —^~|“Oo 

In the ferromagnetic phase, the exponential decay of the typical transverse field allows to define some correlation 
length ^h{0) (see Fig. |l](a)) 


lnh*/P| 


L 


L \9<9c 


oc —- 


L-S-+00 ih{S) 

The growth of the width of the probability distribution of In h (Eq. 1551) 


A 


In/ti — 


(115) 


(ln/iL)2 - (ln/iL)2 oc 


^width 


(116) 
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also allows to define another correlation length ^width- Their divergences near criticality involve the same exponent 
(see Fig. [l](b)) 

t/— 

UdthiO) cx (117) 

0—^9c 

Our numerical results for a = and ct = 1 point towards the same value 

iy{cr = ‘^) ^ 6 ^ = 1) (118) 

in agreement with the finite-size scaling exponent of Eq. IllOl 


G. Comparison with the previous Strong Disorder RG results |lOll 

Let us compare the results described above with the Strong Disorder RG study 0 : 

(i) here the dynamical exponent is found to be z{a) ~ a instead of the value z{a) ~ 1 -|- cr obtained via Strong 

Disorder RG [T^ ■ As a consequence of the relation of Eq. [731 this difference of unity in z is directly related to the 
difference between the magnetic exponent x{a) ~ 0.5 obtained here and x = 1 found in Ref 0 (where more precisely 
the magnetization grows as m{L) oc ). At the level of the renormalization rules, the origin of this difference 

is the following : whereas the present renormalization rules become equivalent to the Strong Disorder RG rules in 
the limit W —>■ -foo (see section riVBI) . they are not equivalent anymore in the critical region. Indeed in the present 
framework, the magnetization is already of order m{L) oc jn the paramagnetic phase, which can be seen as the 

result of Gentral Limit fluctuations. Since the magnetization at criticality cannot be smaller, one obtains the bound 
a; < 1/2. 

(ii) instead of the essential singularity of Eq. | 5 ] 0 , we obtain here conventional power-law scaling laws with finite 
but rather large correlation length exponent ~ 6, and ma gne tic exponent /3 ~ 3. 

In summary, the Strong Disorder RG results obtained in |10l | for the Long-Ranged chain are based on the approxi¬ 
mation that the quantum fluctuations are negligible with respect to disorder fluctuations, whereas the present block 
renormalization for the Dyson version yields that quantum fluctuations are important in the critical region. Further 
work is needed to better understand whether the differences in the results are entirely due to the different approxi¬ 
mations made in the renormalization rules, or whether the Dyson hierarchical version actually changes significantly 
the Long-ranged model and reduces the importance of rare events. 


V. CONGLUSION 

In this paper, we have proposed to study via real-space renormalization the Dyson hierarchical version of the quan¬ 
tum Ising chain with Long-Ranged power-law ferromagnetic couplings J(r) oc and pure or random transverse 

fields. For the pure case, the RG rules have been explicitly solved as a function of the parameter cr, and we have 
compared the critical exponents with previous results of other approaches. For the random case, we have studied 
numerically the RG rules and compared the critical properties with the previous Strong Disorder Renormalization 
approach [lOl| . 

Our conclusion is that the Dyson hierarchical idea initially developed for classical spins is also very useful for 
quantum spin models, in order to derive real-space renormalization rules. In the future, we hope to apply this method 
to other quantum models. 


Appendix A: Reminder on the Mean-Field theory for the pure Long-Ranged model 

1. Thermodynamic exponents 

For the pure quantum Ising chain of Eq. |2l the mean-field theory amounts to look for the uniform magnetization 
m 


< erf > = m 

< tjf > = \/l — 


(Al) 








that minimizes the ground state energy per spin 


(A2) 


easim) = 


Ecs{m) 


— -= -h\/l -m? - 


where Jtot represents the sum of all couplings linked to a given spin 

Jtot = Jl,j 

For instance for the Dyson model of Eq. [181 it is given by 


+00 +00 

jDyson __ \ ^ ji^) _ \ ^ _ _- 

tot / j / j ^_ 2 " 


The minimization of Eq. IA2I 

0 = dmeasim) = h 

yields that the critical transverse field is 


k=0 


\/l — mP 


- Jtotm 


hr - Jf.l 


and that the optimal magnetization reads 

m =0 


for h > hr 


m = W 1 — ( — for h < he 

corresponding to the usual mean-field thermodynamic exponents 

Pmf = 2 
Q-mf = 0 

independently of the short-range or long-ranged nature of the couplings. 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(AS) 


2. Correlation exponents for the short-ranged case 


As recalled in the Introduction, the short-ranged model is characterized by the dynamical exponent (both in mean- 
field and outside mean-field) [l| 

= 1 (A9) 

The correlation length exponent f and the exponent rj of the critical correlation are given by the standard values 

,,SR _ i 

^MF — 2 

vlfF =0 (AlO) 

The upper critical dimension du above which mean-field exponents apply is the value where the mean-field exponents 
satisfy the hyperscaling relation 


leading to 


2 — a = {d + z)v 


du — 


2 — aMF 


.SR 


(All) 


(A12) 


- = 3 
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3. Correlation exponents for the long-ranged case 


For the long-ranged model in dimension d 


-d—0 


J(r) oc r 

the Gaussian fixed-point of the paramagnetic phase h > corresponds to the correlation Q 


(A13) 


/ +°° fldZ ^-too 


dui 


^iik.r+bjt) 


{2t:Y J_^ 27r (/i_/i^)+a;2-h |fc|<" 


(A14) 


where the singular term 1^1“^ replaces the standard Short-Ranged term for cr < 2, with the following consequences : 
(i) the dynamical exponent z is fixed by the anisotropy between and \kY 


^MF — 


(ii) the spatial correlation length exponent v is fixed by the balance between the terms \h — hc\ and |fc|'^ 


LR 

^MF 


1 

cr 


(A15) 


(A16) 


(ii) the temporal correlation length is fixed by the balance between the terms {h — he) and so that the gap 
(corresponding to the inverse temporal correlation length) vanishes as 


Amf = (h - with gMF = 


1 


(A17) 


as in the short-ranged case. By consistency with Eq. lAlSI and IA161 one has of course gMF = 
(iv) the global static susceptibility diverges 

/ -l-oo ^-1-00 ^ 

d‘^rj^ dtC{r,t) = with7MF = l 


(v) at criticality h = he, the correlation reads 

Ceriti[r,t) = / 


ik.', 


^+oo 


du e*‘ 


d<^k = * 


(27r)“ 


27r a;2 -h \kY 


(27r)'’ 


2 | fc|1 


In particular, the spatial correlation at coinciding points follows the power-law 


CeritiY 5 ^ — 


+°° d^k 
-oo (27r)'^2|fc|f 


oc = j.-(.d-2+z+Ti) 


with 


(A18) 


(A19) 


(A20) 


Vmf — 2 — cr 


(A21) 


The upper critical value tJu below which mean-field exponents apply is the value where the mean-field exponents 
satisfy the hyperscaling relation of Eq. lAllI 


2 — cxmf — {d + Zmf)^mf 


(A22) 


leading to 0 






In particular in dimension d = 1, the mean-field region d > d^ corresponds to the region a < with 

2 


(A23) 


(A24) 
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For comparison with the RG procedure described in the text, it is useful to mention the corresponding exponents 


and the magnetic exponent 


’^{cTu = 
z(au = 


1 3 

(7u 2 

^ _ 1 

2 “ 3 


x{cju = 




Pmf _ 1 
= |) 3 


(A25) 


(A26) 


Note that at the upper critical value cr^, one can still use the standard finite-size-scaling valid in the non-mean-held 
region cr > cr„ with the mean-held exponents valid in the mean-held region 0 < tr < cr^. But for 0 < cr < cr^, the 
standard hnite-size-scaling does not hold anymore and is replaced by modihed hnite-size-scaling (see the series of 
recent works (48l - [5ll | and references therein). 
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